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NOMENCLATURE 


a 

inner  radius,  used  as  length  unit  in  this  analysis,  a  =  1 

b 

diameter  ratio  or  normalized  outer  radius 

c 

normalized  depth  of  cracks 

E 

Young’s  modulus  of  the  tube  material,  30x10^  psi 

K 

opening  mode  stress  intensity  factors  (SIF) 

K(Po) ,K(Pi) 

SIF  due  to  uniform  tension  pQ  and  internal  pressure  p^ 

SIF  due  to  crack  face  loading  pc 

Kc(p) ,Kc(pr~2) 

SIF  due  to  pc  =  p,  pc  -  pr“2,  respectively 

Kc( e=0 .9) 

SIF  due  to  pc  corresponding  to  a  90  percent  overstrain 
residual  stress 

N 

number  of  radial  cracks 

r,0 

polar  coordinates  centered  at  the  center  of  the  tube 

rc 

radius  to  the  tip  of  an  0D  crack,  rc  =  b  -  c 

t  * 

wall  thickness  of  the  cylinder,  t  =  b  -  1 

T>T0>Tp 

temperature  at  r,  r  =  a  and  r  =  p,  respectively 

a 

linear  thermal  expansion  coefficient,  6.8x10“^  in./in./°F 

e 

percentage  of  overstrain,  e  =  (p-l)/t,  0  <  e  <  1 

V 

Poisson’s  ratio,  0.3 

P 

radius  of  elastic-plastic  interface  during  pressurization 

ao. 

uniaxial  yield  stress  of  the  tube  material,  170  Ksi 

Q 

Q 

CD 

normal  stress  in  the  radial  and  tangential  direction, 
respectively 
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I  'JTRODUCTION 


An  analytic  method  is  not  available  for  the  computation  of  stress 
intensity  factors  for  multiple-radial  cracks  emanating  from  the  bore  of  a 
thick-wall  cylinder.  The  computation  must  depend  on  various  numerical 
methods.  The  finite  element  method,  using  12-node  quadrilateral, 
isoparametric  elements,  has  been  used  to  obtain  reasonably  accurate  results 
for  a  uniform  array  of  equal  depth  radial  cracks  near  the  bore  of  a  thick-wall 
cylinder  subjected  to  internal  pressure. 1  It  is  common  practice  to 
autofrettage  a  cylinder  so  a  compressive  tangential  stress  is  introduced  near 
the  bore.  This  compressive  residual  stress  reduces  the  tensile  stress  near 
the  bore  caused  by  bore  pressure.  It  increases  the  maximum  internal  pressure 
a  cylinder  can  contain  elastically,  and  extends  the  fatigue  life.  The 
residual  stress  causes  added  difficulties  in  computing  stress  intensity 
factors.  The  thermal  simulation  method2  is  used  to  simulate  the  autofrettage 
residual  stress  by  an  active  thermal  load.  The  combination  of  weight  function 
method  and  the  finite  element  method  developed  in  reference  3  is  used  to 
facilitate  the  numerical  computations  for  various  degrees  of  autofrettage. 


1pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  a  Circular  Ring 
With  Uniform  Array  of  Radial  Cracks  Using  Cubic  Isoparametric  Singular 
Elements,”  Fracture  Mechanics,  ASTM  STP  677,  1979,  pp.  685-699. 

Hussain,  M.  A.,  Pu,  S.  L. ,  Vasilakis,  J.  D. ,  and  O'Hara,  P.,  "Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,"  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 

pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in 
a  Partially  Autof rettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 
National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles,  1981  and 
published  as  a  Technical  Report  ARLCB-TR-81027 ,  1981. 
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While  the  autofrettage  process  produces  favorable  compressive  residual 
stress  near  the  bore,  it  also  yields  a  tensile  residual  stress  near  the  outer 
cylindrical  surface.  This  tensile  stress,  in  addition  to  the  tensile  stress 
produced  by  a  bore  pressure,  causes  a  high  level  of  stress  that  may  enhance 
crack  initiation  or  propagation  from  the  outer  surface  of  the  cylinder.  The 
problem  of  an  externally  cracked,  thick-wall  cylinder  subjected  to  an  internal 
pressure  has  been  studied  by  Kapp^  using  the  finite  element  method  and  by 
racy^  using  the  modified  mapping  collocation  method.  The  K  solution  for  a 
partially  autof rettaged  cylinder  with  external  cracks  was  obtained  recently  by 
Parker.6  The  purpose  of  this  report  is  to  extend  the  method  of  functional 
stress  intensity  factor  recently  developed  in  reference  3  for  internally 

cracked  cylinders  to  externally  cracked,  partially  autof rettaged ,  thick-walled 
cylinders . 

Since  an  application  of  this  study  to  gun  barrels,  some  basic  differences 
are  noted  between  ID  and  OD  cracks  in  cannon  tubes.  A  large  number  of  ID 
cracks  appears  after  only  a  few  rounds  of  firing  due  to  heat  crazing.  The 
number  of  OD  cracks  is  far  less.  The  growth  of  an  ID  crack  tends  to  drive  the 


Pu»  s*  L*  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in 
a  Partially  Autof rettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 
National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles,  1981  and 
published  as  a  Technical  Report  ARLCB-TR-81027 ,  1981. 

Kapp,  J.  A.  and  Eisenstadt,  R. ,  "Crack  Growth  in  Externally  Flawed, 

Autof rettaged  Thick-Walled  Cylinders  and  Rings,"  Fracture  Mechanics,  ASTM  STP 
^677,  1979  ,  pp.  746-756 . 

Tracy,  P.  G. ,  "Elastic  Analysis  of  Radial  Cracks  Emanating  From  the  Outer  and 
Inner  Surfaces  of  a  Circular  Ring,"  Engineering  Fracture  Mechanics.  Vol.  11 
1979  ,  pp.  291-300.  - ~~ - 

aParker ,  A.  P. ,  "Stress  Intensity  and  Fatigue  Crack  Growth  in  Multiply- 
Cracked,  Pressurized,  Partially  Autof rettaged  Thick  Cylinders,"  AMMRC 
TR81-37,  1981. 
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crack  Lip  away  from  the  favorable  compressive  tangential  stress  field  while 
the  growth  of  an  OD  crack  leads  the  crack  tip  inward  into  the  compressive  zone 
of  an  autof rettaged  cylinder.  Therefore,  the  numerical  evaluation  of  K  for 
OD  cracks  will  be  limited  in  this  report  to  a  small  number  of  cracks  with  a 
crack  depth  relatively  small  so  that  the  crack  tip  remains  in  the  tensile 
stress  field. 

FINITE  ELEMENT  METHOD 

In  the  past  decade  various  finite  element  techniques  have  been  developed 
to  successfully  treat  the  singular  stress  field  associated  with  cracks.  The 
use  of  high  order  elements  is  favored  because  of  high  degree  of  accuracy 
obtainable  with  a  small  number  of  elements.  This  prompted  us  to  choose  12- 
node  quadrilateral,  isoparametric  elements  for  the  solution  of  the  problem  of 
multiple-radial  cracks  in  a  thick-wall  cylinder.  A  special  collapsed 
triangular  element/  was  developed  for  linear  elastic  fracture  mechanics 
analysis.  No  sooner  had  the  12-node  quadrilateral,  isoparametric  elements  and 
the  special  crack  tip  elements  been  implemented  in  the  general  purpose  NASTRAN 
program  than  we  found  the  special  computer  program  APES  (Axi symmetric /Planar 
Elastic  Structures),  developed  by  Gifford.^  The  APES  program  was  written 
specifically  for  the  use  of  the  12-node  quadrilateral,  isoparametric  elements. 
It  is  a  user  oriented  program  and  is,  in  fact,  very  easy  to  use.  APES 
consists  of  two  types  of  crack  tip  elements  for  linear  elastic  fracture 

7Pu,  S.  L.  and  Hussain,  M.  A.,  "The  Collapsed  Cubic  Isoparametric  Element  as  a 
Singular  Element  for  Crack  Problems,"  International  Journal  of  Numerical 
Methods  in  Engineering,  Vol.  12,  1978,  pp.  1727-1742. 

^Gifford,  L.  N. ,  Jr.,  "APES-Second  Generation  Two-Dimensional  Fracture 
Mechanics  and  Stress  Analysis  by  Finite  Elements,"  DTNSRDC  Report  4799,  1975. 
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analysis.  The  "enriched"  element  is  deemed  better  than  the  "core"  crack  tip 
element.  We  have  implemented  the  newly  found  collapsed  triangular  elements 
into  APES  as  another  option  for  treating  linear  fracture  structures.  Both 
enriched  elements  and  collapsed  triangular  elements  yield  comparably  accurate 
results  for  stress  intensity  factors.  However,  the  accuracy  loss  due  to 
distortion  in  isoparametric  finite  elements^ > ^  should  be  considered  in 
choosing  particular  crack  tip  elements  in  setting  up  the  finite  element 
idealization  for  a  given  geometrical  configuration.  The  collapsed  elements 
were  used  in  reference  1.  The  accuracy  of  the  results  was  good  in  general 
except  several  values  suffered  accuracy  loss  due  to  excessive  distortion.  The 
enriched  elements  were  used  in  reference  3  for  interior  cracks  in  a  partially 
autof rettaged  cylinder.  Table  I  shows  a  good  agreement  between  the  finite 
element  results* * 3  and  the  trodified  mapping  collocation  results  of  Tracy3  for 
non-autof rettaged  cylinders.  A  similar  agreement  between  results  of  these  two 
methods  for  fully  autof rettaged  cases  is  shown  in  Table  II.  The  irodified 
mapping  collocation  results  in  Table  II  were  obtained  from  data  supplied  by 


XPu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  a  Circular  Ring 

With  Uniform  Array  of  Radial  Cracks  Using  Cubic  Isoparametric  Singular 
Elements,"  Fracture  Mechanics,  ASTM  STP  677,  1979,  pp.  683-699. 

3Pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in 
a  Partially  Autof rettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 
National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles,  1981  and 
published  as  a  Technical  Report  ARLCB-TR-81027 ,  1981. 

^Tracy,  P.  G. ,  "Elastic  Analysis  of  Radial  Cracks  Emanating  From  the  Outer  and 
Inner  Surfaces  of  a  Circular  Ring,"  Engineering  Fracture  Mechanics,  Vol.  11, 
1979,  pp.  291-300. 

^Hopkins,  D.  A.  and  Gifford,  L.  N. ,  "Accuracy  Loss  in  Distorted  Isoparametric 
Finite  Elements,"  DTNSRDC  Structures  Department  Report  M2,  1978. 

^Sickles,  J.  B.  and  Gifford,  L.  N.  ,  "A  Further  Study  of  Accuracy  Loss  in 
Distorted  Isoparametric  Finite  Elements,"  DTNSRDC  Structures  Department 
Report  M-50,  19  79. 
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Parker multiplied  by  the  factor  (2//3)  which  converts  a  result  in  Trescafs 
yield  criterion  to  a  corresponding  result  in  von  Mises’  criterion.  Due  to  the 
success  in  obtaining  accurate  results  for  the  case  of  inner  cracks,  the 
enriched  elements  were  again  used  in  this  report  for  multiple  cracks  emanating 
from  the  outer  surface.  Figure  1(a)  shows  a  typical  finite  element 
idealization.  For  a  very  shallow  crack,  the  slight  modifications  in  finite 
element  meshes  are  shown  in  Figure  1(b). 

TABLE  I.  COMPARISON  OF  iq/pi/irc  FOR  ID  CRACKS  IN  A  CYLINDER  OF  b/a  =  2, 
SUBJECT  TO  BORE  PRESSURE  pf.  FINITE  ELEMENT3  VS.  MODIFIED 
MAPPING  COLLOCATION  (TRACY3). 


1  Crack 

2  Cracks 

4  Cracks 

c 

FE 

MMC 

FE 

MMC 

FE 

MMC 

0. 1 

2.83 

2.80 

2.87 

2.83 

2.83 

2.75 

0.2 

2.83 

2.77 

3.01 

2.96 

2.78 

2.75 

0.3 

2.89 

2.88 

3.28 

3.23 

2.83 

2.77 

0.4 

3.02 

2.99 

3.56 

3.55 

2.95 

2.99 

0.5 

3.19 

3.15 

4.00 

4.00 

3.15 

3.15* 

*This  value  is  based  on  H|  =  1.18  in  a  private  communication  prior 

to  the  publication  of  reference  5.  The  value  =  1.31  in  reference 
3  is  probably  a  typographic  error. 


JPu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in 
a  Partially  Autof rettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 
National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles,  1981  and 
^published  as  a  Technical  Report  ARLCB-TR-81027 ,  1981. 

--Tracy,  P.  G.  ,  "Elastic  Analysis  of  Radial  Cracks  Emanating  From  the  Outer 
and  Inner  Surfaces  of  a  Circular  Ring,"  Engineering  Fracture  Mechanics, 

Vol .  11,  1979,  pp.  291-300. 

1  Parker,  A.  P.  and  Andrasic,  C.  P. ,  "Stress  Intensity  Prediction  For  a 

Multiply-Cracked,  Pressurized  Gun  Tube  With  Residual  and  Thermal  Stresses," 
Army  Symposium  on  Solid  Mechanics  1980  and  supplemental  data  supplied  in 
private  communications,  1981. 
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TABLE  II. 


COMPARISON  OF  iq /-0o/irc  FOR  ID  CRACKS  IN  A  FULLY  AUTOFRETTAGED 
CYLINDER  OF  b/a  =  2.  FINITE  ELEMENT3  VS.  MODIFIED  MAPPING 
COLLOCATION*  (PARKER11). 


4  Cracks 

6  Cracks 

10 

Cracks 

40  Cracks 

c 

FE 

MMC 

FE 

MMC 

FE 

MMC 

FE 

MMC 

0.05 

0.981 

0.967 

0.975 

0.959 

0.957 

0.941 

0.659 

0.670 

0.1 

0.895 

0.871 

0.871 

0.847 

0.817 

0.794 

0.442 

0.432 

0.2 

0.745 

0.722 

0.684 

0.662 

0.575 

0.557 

0.277 

0.271 

0.3 

0.628 

0.609 

0.539 

0.521 

0.417 

0.403 

0.197 

0.191 

-  - -  oaFFi.j.cu  Vjy  raiAd  ill  d  pi  LVdLti  COlIllUUn  LC 3 L  lOtl  LT 

multiplied  by  2/v3  to  convert  from  Tresca  yield  criterion  to  von 
Mises'  criterion. 


THERMAL  SIMULATION 


In  a  partially  autof rettaged  tube,  the  plane-strain  residual  stresses 
based  on  von  Mises'  yield  criterion  for  an  incompressible  material  are  given 

by 13 


°  r  r 

—  {(2  log  -  -  1  +  — ) 

/3  p  b2 


1  <  r  <  p  (1) 


Pr(r) 


✓3 


(P2-Pl)(- 


1 


P  <  r  <  b  (2) 


O 

Pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in 
a  Partially  Autof rettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 
National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles,  1981  and 
published  as  a  Technical  Report  ARLCB-TR-81027 ,  1981. 
iDarker,  A.  P.  and  Andrasic,  C.  P.,  "Stress  Intensity  Prediction  For  a 
Multiply-Cracked,  Pressurized  Gun  Tube  With  Residual  and  Thermal  Stresses," 
Army  Symposium  on  Solid  Mechanics  1980  and  supplemental  data  supplied  in 
private  communications,  1981. 

12Hill,  R.,  The  Mathematical  Theory  of  Plasticity,  Oxford  at  the  Clarendon 
Press,  1950. 
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1  1 


✓3 


{2  log  -  +  1  +  —  -  Pi(—  +  — )  }  1  <  r  <  p 


b2  r2 


(3) 


°e(r)  = < 


/3 


9  1  1 

(p2-Pl)  (—  +  -r) 
b2  r2 


p  <  r  <  b 


where 


b2  p" 

Pi  =  Pi (p)  -  7  (1  -  —  +  2  log  p) 

b2-l  b2 


(4) 


(5) 


In  this  residual  stress  field,  the  presence  of  a  sharp  crack  or  a  notch  may 
cause  geometric  changes  of  the  cylinder  and  redistribution  of  stresses •  The 
stress  redistribution  in  case  of  a  notch  and  the  stress  intensity  factor  in 
case  of  a  sharp  crack  can  not  be  directly  formulated  as  a  boundary  value 
problem  for  finite  element  application.  It  is  verified  in  reference  2  that 
the  following  thermal  stresses  are  present  in  a  thick-wall  cylinder 

Ea(T0-Tp)  r  p2  1  1 

!  -  ((2  log - 1  +  —)  -  Pi(— - -)  }  1  <  r  <  p  (6) 

2(l-v)log  p  P  b2 


b2  r2' 


ar(r)  =  ( 


Ea(T0-Tp) 

(p2  "  Pi) (“r  ~  -r) 


'  2( 1— v) log  p 


b2  r2 


p  <  r  <  b  (7) 


Ea(To-Tp)  r  p2 

-  {2  log  -  +  1  +  — 

2(  l-v)log  p  p  b2 


1  1 
Pl(79  + 


-°e(r)  =  < 


Ea(T0-Tp)  i  i 

; - - -  (P2-Pi)(—  +  — ) 

v  2(  1“ v) log  p  p2  r2 


1  <  r  <  p  (8) 


p  <  r  <  b  (9) 


^Hussain,  M.  A.,  Pu,  S.  L. ,  Vasilakis,  J.  D. ,  and  O’Hara,  P. ,  "Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,”  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 
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if  the  cylinder  is  subjected  to  the  thermal  load 


T(r)  =  J 


-  log  r  1  <  r  <  p 

log  p 


(10) 


T 


P 


p  <  r  <  b 


The  thermal  stresses  and  the  residual  stresses  become  equivalent  if  the 
temperature  gradient  of  the  thermal  load  and  the  material  properties  of  the 
cylinder  satisfy  the  following  relation 


4(  1- v)  oQ 


To  "  Tp  = 


-  log  p 

/3  Ea 


(ID 


where  aQ  is  the  uniaxial  yield  stress  and  T0  or  Tp  may  be  assigned 
arbitrarily. 

The  replacement  of  residual  stresses  by  an  equivalent  thermal  load 
provides  a  method  by  which  the  original  residual  stress  problem  is  transformed 
into  a  boundary  value  problem.  The  application  of  the  thermal  simulation 
method  is  illustrated  in  references  3  and  13.  In  this  report  the  method  of 
thermal  simulation  is  again  used  to  handle  external  radial  cracks  in  a 
residual  stress  field.  Since  r  =  p  divides  the  temperature  into  two  distinct 
fields,  it  is  advisable  to  use  r  *  p  as  a  dividing  arc  in  the  finite  element 
idealization  to  correctly  reflect  the  constant  temperature  field  versus  the 
logarithmically  varying  field. 

JPu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in 
a  Partially  Autof rettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 
National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles,  1981  and 
published  as  a  Technical  Report,  ARLCB-TR-81027 ,  1981. 

'3Pu,  S.  L.  and  Hussain,  M.  A.,  "Residual  Stress  Redistribution  Caused  by 
Notches  and  Cracks  in  a  Partially  Autof rettaged  Tube,"  Journal  of  Pressure 
Vessel  Technology,  Vol.  103,  No.  4,  1981,  pp.  301-306. 
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WEIGHT  FUNCTION  AND  FUNCTIONAL  INTENSITY  FOR  EXTERIOR  CRACKS 


The  finite  element  and  the  thermal  simulation  methods  enable  us  to 
compute  stress  intensity  factors  for  cracks  in  a  cylinder  subjected  to  an 
internal  pressure  and  autofrettage  residual  stresses.  A  combination  of  these 
methods  and  the  weight  function  method  is  developed  in  reference  3.  For  a 
given  crack  geometry,  only  three  finite  element  computations  are  needed  for 
the  three  types  of  loadings  applied  separately.  These  loadings  are  (a)  bore 
pressure  p-^  (not  on  crack  face),  (b)  uniform  tension  pQ  on  OD  and,  (c)  a 
thermal  load  simulating  the  100  percent  overstrain  residual  stresses.  These 
values  of  stress  intensity  factors  are  denoted  by  K(p^),  K(p0),  and  K(e-l). 
Functional  stress  intensity  factors  are  defined  as  Kc ( 1 )  ,  Ke(r”2),  and  K^log 
r)  and  can  be  solved  algebraically  from  K(p^),  K(p G) *  ana  K(e=l).  The  stress 
intensity  factor  for  the  same  crack  geometry  due  to  a  partial  autofrettage 
residual  stress  is  given  by  a  simple  algebraic  equation  in  terms  of  the 
functional  stress  intensity  factors.  The  method  used  for  internal  cracks  in 
reference  3  can  be  equally  applied  to  external  cracks.  The  expressions  for 
external  cracks  given  below  are  slightly  different  from  those  given  in 
reference  3  for  internal  cracks. 

According  to  Bueckner^  and  Rice,  15  a  weight  function  is  a  universal 
function  which  depends  only  on  geometry  and  not  on  loadings.  If  the  stress 


^Pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  For  Radial  Cracks  in 
a  Partially  Autof rettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 
National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles  1981  and 
published  as  a  Technical  Report  ARLCB-TR-8 1027 ,  1981. 
l^Bueckner,  H.  F. ,  "A  Novel  Principle  for  the  Computation  of  Stress  Intensity 
Factors,"  Z.  Agnew  Math.  Mech. ,  Vol.  50,  1970,  pp.  529-546. 
l^Rice,  J.  R.  ,  "Some  Remarks  on  Elastic  Crack-Tip  Stress  Fields,"  Int.  Journal 
of  Solids  and  Structures,  Vol.  8,  1972,  pp.  751-758. 
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intensity  factor  ^  and  displacement  field  ^  associated  with  load  system 
1  are  known,  the  weight  function  for  the  cracked  geometry  is 

.  H  3u(l)(x,y,c) 

h  - - (12) 

2K(D(c) 

where  H  =  E  for  plane  stress  and  H  =  E/(l-v2)  for  plane  strain.  Once  h  is 
determined,  the  stress  intensity  factor  induced  by  any  other  symmetric  load 
system  t  and  f  is  given  by 

K  =  /  fh  dr  +  /  f *h  dA  (13) 

r  a 

where  t  is  the  stress  vector  acting  on  boundary  f  around  the  crack  tip  and  f 
is  the  body  force  in  region  A  defined  by  T.  For  our  radial  cracks  in  a 
cylinder,  f  s  0  is  used  and  equation  (13)  can  be  reduced  to 

H  c  3v 

K  -  -  /  pc(x)  —  dx  (14) 

K*  o  3c 

where  c  is  the  crack  depth  and  x  is  a  distance  both  measured  from  the  outer 
cylindrical  surface  with  x  and  r  related  by 

x  =  b  -  r  (15) 

K*  and  v  are  the  stress  intensity  factor  and  the  crack  opening  displacement 
associated  with  a  certain  loading,  while  K  is  the  stress  intensity  factor  for 
the  applied  symmetric  loading  and  pc(x)  is  the  equivalent  crack  face  loading 
reduced  from  the  applied  loads  by  superposition  of  a  solution  for  a  crack- 
free  body  under  the  same  loads. 

For  a  uniform  tension  p0  on  the  OD  and  an  internal  pressure  (note  pQ 
or  p-L  not  acting  on  crack  faces),  the  equivalent  crack  face  loadings  are 

respectively  j 

i 

! 
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b2Po 

pc(x)  **  ~Z —  [1  +  (b-x)-2]  (16) 

b2-l 

Pi 

Pc(x)  “  [1  +  b2(b-x)-2]  (17) 


The  crack  face  loading  equivalent  to  a  100  percent  autofrettage  residual 
stress  is 

pc(x)  =*  —  [{2  -  Pi(b)}  -  Pi(b)(b-x)-2  +  2  log(b-x)  ]  (18) 

/3 


Substituting  into  equation  (14)  from  equations  (16),  (17),  and  (18)  we  have 

(19) 


K(p0)  b2  Kc(l)  b2  Kc(r"2) 

+ 


r\  — O 

Po/tTC*  k  1  /ttc  ^  “1  /tic 
K(Pi)  1  Kc(l)  b2  Kc(r"2) 

Pi/i^ "  b2-l  /i£  b2-l 


(20) 


Kc(e-l) 


[{2  -  P1(b)}Kc(l)  -  Pi(b)Kc(r-2)  +  2Kc(log  r) ]  (21) 


00/TTC 


'3ttc 


where  the  functional  intensity  factors  are  defined  by 


H  c  3v 

Kc(D  -  "  /  “  dx 

K*  o  3c 


(22) 


9  H  .c  .  3v 

Kc(r-2)  =  -  /  (b-x)-2  —  dx 
K*  o  3c 


(23) 


H  c  3v 

Kcdog  r )  -  -  /  [  log(b-x)  ]  —  dx 

K*  o  3c 


(24) 


For  a  given  crack  configuration,  we  can  solve  for  these  functional 
intensity  factors  from  equations  (19)  to  (21)  with  the  left  hand  sides  of 
these  equations  known  from  finite  element  computations*  Once  the  functional 
intensity  factors  are  known,  it  is  a  simple  matter  to  compute  stress  intensity 
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factors  for  the  same  crack  due  to  residual  stresses  corresponding  to  a  partial 
autof rettage ,  provided  that  OD  cracks  are  in  the  region  beyond  r  =  p.  The 
tangential  residual  stress  in  this  region  is 

ao 

oq(x)  =  —  [p2  -  Pi(p)][b“^  +  (b-x)“2]  ,  0  <  x  <  (l-e)t  (25) 

/3 


Using  Eq.  (25)  as  the  crack  face  loading  in  equation  (14),  we  have 

Kc(e)  i  p2  -  Pi(p)  Kc(l)  o  Kc(r“2) 

- =  —  [ - - - ]  [ -  +  b2  - ]  ,  c  <  (l-e)t 

0o/ttc  /3  b  /ttc 

Comparing  equation  (26)  with  equation  (20)  we  can  write 

Kc(e)  i  K(pi) 

- - —  =  —  (p2  -  1  -  2  log  p)  - 

o0/ttc  /3  Pi/lTC 


(26) 


(27) 


It  is  clear  that  the  stress  intensity  factor  due  to  partial  autofrettage 
residual  stress  is  a  multiplication  factor,  depending  on  degree  of 
autofrettage,  times  the  stress  intensity  factor  due  to  internal  pressure. 
This  was  also  observed  by  Parker. ^ 

In  case  the  crack  tip  depth  c  is  larger  than  (l-e)t,  the  crack  face 
loading  must  be  expressed  in  the  form 

/  °o 

—  [p2  ~  Pi(p)][b-2  +  (b-x)-2]  0  <  x  <  (l-e)t  (28a) 

/3 


PC(X>  = 

°o  9 

—  [{2  -  Pi(p)}  -  Pi(p)(b-x)“2  +  2  log(b-x) ] 

'  /3 

(  1-e) t  <  x  <  t  (28b) 


^Parker,  A.  P.  and  Andrasic,  C,  P,  ,  "Stress  Intensity  Prediction  For  a 
Multiply-Cracked,  Pressurized  Gun  Tube  With  Residual  and  Thermal  Stresses, 
Army  Symposium  on  Solid  Mechanics  1980  and  supplemental  data  supplied  in 
private  communications,  1981. 
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n  this  case  neither  equation  (26)  nor  (27)  is  valid.  However,  if  c/t  =  1  — e 
+  6  and  6  is  small,  we  may  use  equation  (26)  or  (27)  to  obtain  an  approximate 
stress  intensity  factor.  The  final  stress  intensity  factor  is  the  sum  of  the 
approximate  and  a  correction  Kg.  The  correction  stress  intensity  factor  Kg  is 
calculated  from 


1  H 
/3 tc  K* 


X  l-e+S)  t 

(l-e)t 


Pc(x) 


dv 

—  dx 
9c 


(29) 


where  pc(x)  is  the  difference  between  equations  (28b)  and  (28a)  or 

pc(x)  =  (1-2  log  p)  -  p2(b-x)-2  +  2  log(b-x)  (30) 

and  v  is  an  approximation  based  on  the  crack  face  displacement  of  Westergaard 
near  field  solution 

2K*  25  1 /2 

v(5)  -  —  (— )  (31) 

H  IT 

where  5  is  a  distance  measured  from  the  crack  tip  outward  radially. 

Carrying  out  the  integration  of  equation  (29)  using  equations  (30)  and  (31)  we 
have 

1  2 

- —  =  •— :  -  {(1-2  log  p)(I1+I1')  -  p2(I2+I2')  +  2(13+13')}  (32) 

a0/7Tc  /  3ttc  ^ 


where 


2 

Ii  =  2 /fit  ,  Ii»  -  —  (6t)3/2 

3c 

1  ^6t  ,  ,  _ _ _ 

12  -  —  [ -  +  rc  tan  1  J&tT rc ] 

rc  p 

l  /  6t 

^2?  =  “  [-  -  +  rc“l/2  tan“  1/6"t/r>  ] 

c  p 


(33) 

(34) 


(35) 
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I3  =  2[(-2+log  p)/SF  +  2/r^T  tan”  V <St/ rc ] 


(36) 


2  2 

I3'  =  —  [ ( 6t) 3/2iog  p  +  -  (3rc-6t)/6t  -  2rc^/2  tan-1/6t7r^]  (37) 

Jc  3 

with  the  abbreviation  rc  =  b  -  c. 

NUMERICAL  RESULTS  FOR  THE  CASE  OF  b  =  2,  a  *  1 

Before  presenting  results  for  external  multiple  cracks,  several  graphs 
supplement  to  reference  3  for  interior  cracks  are  added  here.  Figure  2  is 
stress  intensity  factor  due  to  residual  stresses  in  a  fully  autof rettaged 
cylinder  as  a  function  of  crack  depths  for  several  values  of  N  from  2  to  40. 
Stress  intensity  factors  are  negative  due  to  compressive  crack  face  loadings. 
Similarly,  the  ordinates  of  Figures  5  and  6  of  reference  3  should  be  negative. 
Figure  3  shows  the  decrease  in  stress  intensity  factor  due  to  an  increased 
degree  of  autofrettage  for  a  single  internal  crack. 

The  stress  intensity  factor  due  to  combined  internal  pressure  and 
autofrettage  residual  stresses,  Figure  4,  remains  highest  for  two 
diametrically  opposite  inner  cracks  and  decreases  raonotonically  as  the  number 
of  cracks  increases  for  various  crack  depths  and  various  degrees  of 
autofrettage . 

For  external  multiple  cracks,  the  dimensionless  stress  intensity  factors 
are  shown  as  a  function  of  crack  depths  in  Figure  5  for  a  non— autof rettaged 
cylinder  and  in  Figure  6  for  a  fully  autof rettaged  cylinder  for  several  values 

o 

Pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in 
a  Partially  Autof rettaged  Thick~Wall  Cylinder,  Presented  at  the  Fourteenth 
National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles,  1981  and 
published  as  a  Technical  Report  ARLCB-TR-81027 ,  1981. 
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of  N  from  1  to  20.  A  comparison  of  our  finite  element  results  with  Tracy’s 
modified  mapping  collocation  results^  is  given  in  Table  III  for  cracks 
emanating  from  the.  outer  surface  of  a  non-autof rettaged  cylinder  subjected  to 
an  internal  pressure.  The  agreement  is  not  as  good  as  shown  in  Table  I  for 
cracks  near  the  inner  surface.  Previous  finite  element  results  were,  in 
general,  higher  than  MMC  results  in  Table  I,  but  the  opposite  appears  in  Table 
III.  The  finite  element  results  using  8-node  quadratic  quadrilateral  elements 
with  collapsed  triangles  as  singular  crack  tip  elements  for  one  external  crack 
in  a  non-autof rettaged  cylinder  subjected  to  internal  pressure  given  by  Kapp^ 
are  also  shown  in  Table  III.  Better  agreement  is  seen  between  references  4 
and  5.  This  suggests  an  idealization  with  more  refined  elements  than  used  in 
Figure  1  may  be  needed  for  better  accuracy  for  external  cracks.  To 
demonstrate  the  method  it  is  deemed  enough  to  use  the  coarse  meshes  of  Figure 
1.  For  partially  autof rettaged  cylinders  Parker’s  results* **  are  based  on 
results  of  non-autof rettaged  cylinders. ^  The  agreement  between  finite  element 
and  modified  mapping  collocation  for  an  autof rettaged  case  can  only  be 
expected  to  be  as  good  as  that  for  the  non-autof rettaged  case. 


^Kapp,  J.  A.  and  Eisenstadt,  R. ,  ’’Crack  Growth  in  Externally  Flawed, 

Autof rettaged  Thick-Walled  Cylinders  and  Rings,”  Fracture  Mechanics,  ASTM 
STP  677,  1979,  pp.  746-756. 

^Tracy,  P.  G.  ,  ’’Elastic  Analysis  of  Radial  Cracks  Emanating  From  the  Outer 
and  Inner  Surfaces  of  a  Circular  Ring,”  Engineering  Fracture  Mechanics,  Vol. 
11,  1979,  pp.  291-300. 

**Parker,  A.  P.  and  Andrasic,  C.  P. ,  ’’Stress  Intensity  Prediction  for  a 

Multiply-Cracked,  Pressurized  Gun  Tube  With  Residual  and  Thermal  Stresses,” 
Army  Symposium  on  Solid  Mechanics  1980  and  supplemental  data  supplied  in 
private  communications,  1981. 
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TABLE  III.  COMPARISON  OF  iq/pi/irc  FOR  OD  CRACKS  IN  A  CYLINDER  OF  b/a  =  2 
SUBJECTED  TO  A  BORE  PRESSURE  pi.  FINITE  ELEMENT  (THIS  STUDY) 
VS.  MODIFIED  MAPPING  COLLOCATION  (TRACY5). 


1  Crack 

2  Cracks 

3  Cracks 

4  Cracks 

c 

FE 

FE4 

MMC 

FE 

MMC 

FE 

MMC 

FE 

MMC 

0. 1 

0.756 

0.806 

0.800 

0.758 

0.813 

0.753 

0.807 

0.747 

0.807 

0.2 

0.838 

0.917 

0.913 

0.846 

0.913 

0.830 

0.880 

0.802 

0.867 

0.3 

0.960 

1.056 

1.04 

0.977 

1.06 

0.922 

1.00 

0.874 

0.947 

Functional  stress  intensity  factors  are  obtained  using  equations  (19), 
(20),  and  (21),  and  the  results  are  plotted  in  Figures  7,  8,  and  9.  These 
results  are  the  basis  for  the  computation  of  stress  intensity  factors  for 
multiple  exterior  cracks  for  various  degrees  of  autof rettage .  For  example, 
given  N  =  4,  readings  taken  from  Figures  7  and  8  are  Kc(p)/(p/iic)  =  1.12  and 
Kc(  pr"2) /( p/rTc)  =  0.320  for  c/t  =  0.2  and  K^(  p) /( p/ric)  *  1.18,  Kc(  pr~2) /( p/ric) 
-  0.360  for  c/ t  =  0.3.  For  e  =  0.8,  the  stress  intensity  factor  due  to  the 
autofrettage  residual  stress  can  be  computed  from  equation  (26)  which  gives 
K(e=0.8) /(Oq/tTc)  =  0.492  for  c/t  =  0.2.  If  equation  (27)  is  used, 

K( Pi) /( pi/irc)  and  the  multiplication  factors  are  0.802  and  0.615  respectively, 
then  K(e=0.8) /(Oq/tTc)  =  0.493.  For  c/t  -  0.3,  the  crack  tip  enters  the  region 
r  <  1.8  (inside  the  elastic-plastic  radius).  Equation  (26)  gives  an 

^Tracy,  P.  G. ,  "Elastic  Analysis  of  Radial  Cracks  Emanating  From  the  Outer  and 
Inner  Surfaces  of  a  Circular  Ring,"  Engineering  Fracture  Mechanics,  Vol.  11, 
1979,  pp.  291-300. 
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approximate  value,  K(e-0.8) /(o0/tk*)  =  0.537.  Using  equation  (32)  we  compute 
Ks/Coo/Tc)  =  -0.05.  The  final  result  for  c/t  =  0.3  is  Kc(e=0.8) /(o0/t\c)  = 
0.487.  A  finite  element  computation  for  the  case  e  =  0.8,  c/t  =  0.3  was 
performed  and  the  result  is  K/(a0/rrc)  =  0.486. 

The  dimensionless  stress  intensity  factor  as  a  function  of  degree  of 
autofrettage  is  shown  in  Figure  10  for  selected  values  of  N  and  c/t.  In  Figure 
11  the  same  is  shown  as  a  function  of  the  c/t  for  the  single  crack  case  for 
several  values  of  e.  The  stress  intensity  factor  drops  considerably  from  the 
fully  autofrettaged  condition  to  a  60  percent  or  70  percent  autofrettage  for 
small  values  of  c/t.  In  a  similar  graph  for  a  shallow  interior  crack,  Figure 
3,  the  stress  intensity  factor  increases  slightly  from  the  lowest  value  corre¬ 
sponding  to  the  fully  autofrettaged  condition  to  60  or  70  percent  autofrettage. 
This  suggests  that  a  partially  autofrettaged  cylinder  may  sometimes  have  the 
advantage  over  a  fully  autofrettaged  one. 

For  a  combination  of  residual  stresses  and  internal  pressure,  the  stress 
intensity  factor  is  simply  an  algebraic  sum.  The  stress  intensity  factor 
normalized  by  Kq  =  a0(r=b)/ifc  is  shown  as  a  function  of  c/t  in  Figure  12  for 
various  combinations  of  residual  stress  and  internal  pressure  p^  =  0o/f  for  the 
cases  of  two  and  twenty  cracks.  The  stress  intensity  factor  normalized  by 
°0/^  ls  g^en  as  a  function  of  N  in  Figure  13.  It  is  clear  that  the  stress 
.ntensity  factor  is  largest  for  N  =  2  under  various  combinations  of  residual 
stresses  and  internal  pressures.  The  stress  intensity  factor  is  monotonically 
decreasing  as  the  number  of  cracks  increases  from  N  =  2.  For  shallow  cracks, 
i.e.  c/t  <  0.1,  the  interaction  among  cracks  becomes  weak.  The  stress 
intensity  factor  for  N  =  2  is  nearly  the  same  as  for  N  =  1. 
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APPROXIMATIONS  FOR  SHALLOW  CRACKS 

The  finite  element  method  has  some  difficulty  in  calculating  stress 
intensity  factors  for  very  shallow  ID  or  OD  cracks  due  to  the  small  size  of 
crack  tip  elements  required  to  solve  the  shallow  crack  problem.  An  approximate 
approach  is  to  extend  the  graphs  of  stress  intensity  factor  versus  crack  depth 
c  to  small  values  of  c  using  the  knowledge  that  the  limiting  value  of 
MPc)/^  approaches  1.12  oq  as  c  ♦  0  where  oq  =  aQ(r=a)  for  inner  cracks  and 
°0  =  ae(r=b)  for  external  cracks.  The  dotted  lines  in  Figures  7,  8,  and  9 
epresent  the  extensions.  The  curves  in  the  region  0  <  c/t  <  0.05  in  Figures 
2,  3,  and  4  in  reference  3  are  also  such  extensions.  Taking  readings  of 
functional  stress  intensity  from  those  graphs  and  using  the  method  described  in 
unis  report  and  in  reference  3,  we  can  compute  stress  intensity  factors  for 
various  values  of  c,  including  small  values,  for  various  values  of  N. 

Another  approximate  method  for  shallow  cracks  suggested  by  Underwood  and 
Throoplfi  assumes  the  crack  face  loading  be  linear  varying  and  uses  the  stress 
intensity  solutions  of  Bueckner8 * * * * * * * * 17  for  uniform  pressure  and  Tada18  and  Benthera19 


8Pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in 
a  Partially  AutoTrettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 

National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles,  1981  and 
published  as  a  Technical  Report  ARLCB-TR-81027 ,  1981. 

Underwood,  J.  H.  and  Throop,  J.  F.,  "Surface  Crack  K-Estimates  and  Fatigue 

Life  Prediction,"  ASTM  STP  687,  1979,  pp.  195-210. 

Bueckner,  H.  F.,  "Some  Stress  Singularities  and  Their  Computation  by  Means  of 

Integral  Equations,"  Boundary  Problems  in  Differential  Equations.  Ed.,  R.  E. 

Langer,  University  of  Wisconsin  Press,  Madison,  WI,  I960. 

*»Tada,  H.,  Paris,  P.  C.  and  Irwin,  G.  R. ,  The  Stress  Analysis  of  Cracks 

Handbook,  Del  Research  Corporation,  Heller town,  PA,  1973. 

Benthera,  J.  P.  and  Koiter,  W.  T.,  "Asymptotic  Approximations  to  Crack 

Problems,  Methods  of  Analysis  of  Crack  Problems,  Ed.,  G.  C.  Sih,  Noordhoff 

International  Publishing,  1972. 
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for  linear  varying  pressure.  For  a  single  crack  at  inner  surface  the  formula 
is  given  by 

K x  =  [1.12  oe(r=a)  -  0.68  (oe(r=a)  -  oe(r=rc)  }]  /nc  (38) 

A  similar  equation  for  a  single  crack  at  outer  surface  is 

Kj  =  [1.12  oe(r=b)  -  0.68  (oe(r=b)  -  oe(r=rc)  }] /ttc  (39) 

These  equations  give  good  approximations  for  more  than  one  crack  for  small 
values  of  c. 

We  have  worked  out  several  examples  using  both  the  functional  intensity 
factor  approach  and  the  linear  varying  stress  approach  and  compared  the  results 
with  available  results  by  other  methods.  The  comparison  is  given  in  Table  IV 
for  inner  cracks  and  Table  V  for  outer  cracks.  The  agreement  between  the  two 
approximate  approaches  is  good  for  small  c/t  (N  may  be  large).  If  c/t  is  not 
very  small  and  if  N  is  large  then  the  functional  intensity  approach  is  better. 
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TABLE  IV.  COMPARISON  OF  K/Kq*  FOR  SHALLOW  ID  CRACKS  IN  A  CYLINDER  OF  b/a  =  2 
SUBJECTED  TO  EITHER  A  BORE  PRESSURE  OR  A  RESIDUAL  STRESS 
CORRESPONDING  TO  e  3  1  OR  e  =  0.6. 


• 

All  N 

N  -  4 

N  =  20 

N  =  18 

c/ 1 

Loading 

Eq.  (38) 

FE3 * * * * * * 

MMC11 

FE3 

MMC11 

Pi 

1.113 

1.109 

1.  106 

1.084 

1.101 

0.01 

e  =  1 

1.075 

1.067 

1.068 

1.046 

1.063 

e  =  0.6 

0.927 

0.919 

0.902 

Pi 

1.094 

1.077 

1.076 

0.980 

1.010 

0.04 

e  =  1 

1.011 

0.992 

0.994 

0.907 

0.930 

e  =  0.6 

0.866 

0.848 

0.777 

Pi 

1.061 

1.058 

1.035 

0.778 

0.796 

0. 10 

e  =  1 

0.896 

0.892 

0.871 

0.642 

0.657 

e  =  0.6 

0.756 

0.752 

0.539 

*Kq  =  -  pi  /tt"c  for  pressure  pi  acting  on  bore  and  crack  face, 
Kq  =  a o/ttc”  for  autofrettage  residual  stress. 


3Pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in 

a  Partially  Autof rettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 

National  Symposium  on  Fracture  Mechanics,  UCLA,  Los  Angeles,  1981  and 

published  as  a  Technical  Report  ARLCB-TR-81027 ,  1981. 

^-Parker,  A.  P.  and  Andrasic,  C.  P. ,  "Stress  Intensity  Prediction  for  a 

Multiply-Cracked,  Pressurized  Gun  Tube  With  Residual  and  Thermal  Stresses," 

Array  Symposium  on  Solid  Mechanics  1980  and  supplemental  data  supplied  in 

private  communications,  1981. 
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TABLE  V.  COMPARISON  OF  K/K Q*  FOR  SHALLOW  OD  CRACKS  IN  A  CYLINDER  OF  b/a  =  2 
SUBJECTED  TO  EITHER  A  BORE  PRESSURE  p-j^  OR  A  RESIDUAL  STRESS 
CORRESPONDING  TO  e  =  1  OR  e  -  0.6. 


c/t 

Loading 

Eq.  (39) 
All  N 

FE  (This 
N  =  1 

Study) 

N  *  20 

FEa 

N  =  1 

Pi 

1.123 

1.113 

1.102 

1. 124 

0.01 

e  =  1 

0.69 

0.68 

0.67 

e  =  0.6 

0.265 

0.263 

Pi 

1.134 

1.104 

1.052 

1.132 

0.04 

e  =  1 

0.672 

0.654 

0.608 

e  =  0.6 

0.271 

0.263 

0.251 

Pi 

1.157 

1.  128  A 

0.944 

1.208 

0.1 

e  =  1 

0.636 

0.629 

0.525 

e  =  0.6 

0.276 

0.269 

0.225 

*Ko  -  “  Pi/irc  for  bore  pressure  Pi  and  Kq  =  a0/irc  for  autofrettage 


residual  stress. 

AMMC  result  for  this  case  is  K/Kq  =  1.20. 5 


AKapp,  J.  A.  and  Eisenstadt,  R. ,  "Crack  Growth  in  Externally  Flawed, 

Autof rettaged  Thick-Walled  Cylinders  and  Rings,"  Fracture  Mechanics.  ASTM  STP 
,.677 ,  1979,  pp.  746-756.  - 

Iracy,  P.  G. ,  Elastic  Analysis  of  Radial  Cracks  Emanating  From  the  Outer  and 
Inner  Surfaces  of  a  Circular  Ring,"  Engineering  Fracture  Mechanics.  Vol.  11 
1979,  pp.  291-300. 
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CONCLUSIONS 


The  method  developed  in  a  previous  paper3  for  multiple- radial  cracks  at  the 
inner  surface  of  an  autof rettaged  thick- wall  cylinder  can  be  applied  to 
multiple  cracks  at  the  outer  surface  of  the  cylinder  with  slight  modification. 

A  simple  modification  factor  depends  on  the  degree  of  autof rettage,  and  can  be 
used  to  obtain  stress  intensity  factors  due  to  partially  autof rettaged  residual 
stress  and  internal  pressure,  provided  the  OD  cracks  are  not  beyond  the  elastic 
and  plastic  interface  (during  the  autof rettage  pressurization) . 

The  agreement  between  numerical  results  from  finite  element  and  modified 
mapping  collocation  for  the  exterior  crack  case  is  not  as  good  as  for  interior 
cracks.  In  virtue  of  the  results  in  reference  4,  the  accuracy  of  the  finite 
element  results  using  12-node  cubic  elements  probably  can  be  improved  with  an 
idealization  consisting  of  more  elements  than  those  used  in  this  study. 

As  in  the  case  of  interior  cracks,  the  cylinder  with  two  diametrically 
opposed  cracks  is,  in  general,  the  weakest  configuration.  For  more  than  two 
cracks,  the  stress  intensity  factor  decreases  as  the  number  of  external  cracks 
increases.  The  autof rettage  process  will  prolong  the  useful  life  of  a  pressure 
vessel  with  radial  cracks  at  the  inner  surface,  but  will  have  an  adverse  effect 
on  radial  cracks  at  the  outer  surface.  If  100  percent  autof rettage  causes  the 
concern  of  fracture  failure  shifted  from  inner  cracks  to  outer  cracks,  then  the 


3Pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  for  Radial  Cracks  in  a 
Partially  Autof rettaged  Thick-Wall  Cylinder,"  Presented  at  the  Fourteenth 
National  Symposium  on  Fracture  Mechanics ,:  UCLA,  Los  Angeles,  1981  and 
published  as  a  Technical  Report  ARLCB-TR-81027 ,  1981. 

HKapp,  J.  A.  and  Eisenstadt,  R. ,  "Crack  Growth  in  Externally  Flawed, 

Autof rettaged  Thick-Walled  Cylinders  and  Rings,"  Fracture  Mechanics,  ASTM  STP 
677,  1979,  p.  746-756.  - 
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degree  of  autofrettage  should  be  reduced.  A  slight  decrease  in  the  degree  of 
autofrettage  will  slightly  increase  the  stress  intensity  factor  at  interior 
cracks,  and  at  the  same  time  will  considerably  decrease  the  stress  intensity 
factor  at  exterior  cracks.  It  is  desirable  to  obtain  the  optimal  degree  of 
autofrettage  to  best  suit  the  crack  configuration  and  loading. 
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CRACK  TIP 


(a) 


Figure  1(a).  A  typical  finite  element  idealization. 


Figure  1(b).  Idealization  for  very  shallow  cracks. 
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C/t 
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Kc  (C- 1) 

o0/nc 


Figure  2.  Stress  intensity  factor  as  a  function  of  c/t  for  N  ID 
cracks  in  a  fully  autofrettaged  cylinder  of  b/a  =  2. 
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Figure  3.  K/K0  as  a  function  of  c/t  for  a  single  crack  at  inner 

surface  of  a  cylinder  of  b/a  =  2  subjected  to  an  internal 
pressure  pi  =  a0/f,  f  =  3,  and  residual  stresses 
corresponding  to  various  degrees  of  autofrettage ,  e. 
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30 


Figure  4.  K/o0Ac"  for  N  radial  cracks  at  inner  surface  of  a  cylinder  of  b/a  =  .2 

subjected  to  combined  internal  pressure  pj  =  cr0/f  and  residual  stresses 
corresponding  to  given  degrees  of  autofrettage ,  e. 


K 


C/t 

Figure  5.  K/K0  as  a  function  of  c/t  for  N  radial  cracks  at  outer  surface 
of  a  cylinder  of  b/a  =  2  subjected  to  an  internal  pressure  p^. 
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Kc(  C- 1 ) 

<Vnc 


Figure  6.  YL/oqt/vz  as  a  function  of  c/t  for  N  radial  cracks  at  outer 
surface  of  a  fully  autofrettaged  cylinder  of  b/a  =  2. 
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Figure  7. 


b=2 


j - 1 - 1 - 1 _ i _ i 

0.1  0.2  0.3 

c/t 

K?(p)/pAc’  as  a  function  of  c/t  for  N  external  radial  cracks 
with  constant  crack  face  loading. 
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Kp(PclX)) 


P/n<3 


Figure  8.  Kc/p/ttc  vs.  c/t  for  N  external  radial  cracks  with  crack  face 
loading  pc (x)  =  pfb-x)'2. 
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Figure  9.  Kc/p/iT(;  vs.  c/t  for  N  external  radial  cracks  with  crack  face 
loading  pc(x)  =  plog(b-x). 
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Kc  ( £) 


e 


Figure  10.  Kc(e)/a0/iTc  as  a  function  of  e  for  a  cylinder  of  b/a  =  2  for 
several  selected  values  of  N  and  c  (c  =  c/t  since  t  =  1) . 
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K 


C/t 


Figure  11.  K/KQ  as  a  function  of  c/t  for  a  single  externally  cracked 

cylinder  of  b/a  =  2  subjected  to  combined  internal  pressure 
Pi  =  crQ/f,  f  =  3,  and  residual  stresses  corresponding  to  a 
degree  of  autofrettage  e  given. 
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Figure  12.  K/K0  as  a  function  of  c/t  for  N  externally  cracked  cylinders 
of  b/a  =  2  subjected  to  several  combinations  of  internal 
pressure  pi  =  o0/f  and  residual  stresses. 
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Figure  13.  K/Gq/ttc  for  N  radial  cracks  at  outer  surface  of  a  cylinder  of 
b/a  =  2  subjected  to  combined  internal  pressure  pt  =  a0/f, 
where  f  =  1.5  except  otherwise  indicated,  and  residual 
stresses  corresponding  to  given  degrees  of  autofrettage  e . 
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